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Abstract. Metallicities ([Fe/H]) from low resolution spectroscopy obtained with the Very Large Telescope (VLT) 
are presented for 98 RR Lyrae and 3 short period Cepheids in the bar of the Large Magellanic Cloud. Our metal 
abundances have typical errors of ±0.17 dex. The average metallicity of the RR Lyrae stars is [Fe/H] = — 1.48 ± 
0.03±0.06 on the scale of Harris (1996). The star-to-star scatter (0.29 dex) is larger than the observational errors, 
indicating a real spread in metal abundances. The derived metallicities cover the range —2.12 <[Fe/H]< —0.27, 
but there are only a few stars having [Fe/H]> —1. For the ab— type variables we compared our spectroscopic 
abundances with those obtained from the Fourier decomposition of the light curves. We find good agreement 
between the two techniques, once the systematic offset of 0.2 dex between the metallicity scales used in the 
two methods is taken into account. The spectroscopic metallicities were combined with the dereddened apparent 
magnitudes of the variables to derive the slope of the luminosity-metallicity relation for the LMC RR Lyrae stars: 
the resulting value is 0.214± 0.047 mag/dex. Finally, the 3 short period Cepheids have [Fe/H] values in the range 
—2.0 <[Fe/H]< —1.5. They are more metal-poor than typical LMC RR Lyrae stars, thus they are more likely 
to be Anomalous Cepheids rather than the short period Classical Cepheids that are being found in a number of 
dwarf Irregular galaxies. 
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• '-h . 1. Introduction tion off the Zero Age Horizontal Branch (ZAHB; Sandage 

/S ; 1990). 

^ ; The Large Magellanic Cloud (LMC) is a first step in The precise form of the M V (RR)- [Fe/H] relation- 
our knowledge of the extragalactic astronomical dis- smp i s still matter of debate; and while growing con- 
tance scale. A variety of distance estimates to the LMC sensus is being reached on its zero point ((3 ~ 0.6 mag- 
have been published in recent years (for references, at [Fe/H]=-1.5, Cacciari & Clementini 2003, and refer- 
see Clementini et al. 2003a). Among those based on ences therein), there is still disagreement on the value of 
Population II indicators, most are tied to the value of the s i ope since literature values range from 0.30 to 0.18- 
the absolute visual magnitude of the RR Lyrae variables, .20 mag/dex (Cacciari 1999, Carretta et al. 2000). It has 
that is the radial pulsators within the classic instability even been SU g ge sted that a may not be unique on the 
strip on the core He-burning phase of the metal-poor stars metallicity range spanned by the Galactic globular clus- 
(the so-called horizontal branch, HB). The mean intrinsic t ers (GGCs) since two different linear relationships would 
luminosity of the RR Lyrae stars depends on metallic- be followed by the metal-poor ([Fe/H]< -1.5) and the 
ity according to a relation often assumed to be linear: me t a l-rich ([Fe/H]> -1.5) Galactic and LMC cluster RR 
M V (RR) = a[Fe/H] + /3 (Sandage 1981a,b), and, at fixed Lyrae gtars ( Rey et a i. 2 000, Caputo et al. 2000). 
metal abundance, it is also slightly dependent on evolu- Sandage (1993a) found a rather steep slope of the 

My (RR)— [Fe/H] relationship a ~ 0.30 mag/dex from the 



Send offprint requests to: R.G. Gratton from the so-called period-shift analysis of the RR Lyrae 

* Based on data collected at the European Southern stars in a number of Galactic globular clusters. 
Observatory, proposal numbers 62.N-0802, 66.A-0485, and The theoretical evolutionary and pulsation models of 

68.D-0466 horizontal branch (HB) stars favour instead a shallower 
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slope a ~ 0.2 mag/dex, and also provide evidence for 
a non-linearity that makes the relation steeper for more 
metal rich stars, as recently reviewed and summarized by 
Cacciari & Clementini (2003). A mild slope was also found 
from the Baade-Wesselink method applied to Galactic 
field RR Lyrae stars (Fernley et al. 1998), and from the 
determination of the magnitude of the HB in globular clus- 
ters of M31 (Rich et al. 2001), However both these studies 
do not find clear evidence for a break at [Fc/H]=— 1.5, 
thus raising concerns on the actual universality of the HB 
luminosity-metallicity relationship. 

The slope of the luminosity-metallicity relation of the 
RR Lyrae stars can provide clues on the Galactic globular 
clusters relative ages and on the time scale of the Galactic 
halo formation and early evolution. In fact, given the al- 
most constant difference in luminosity between Turn-Off 
and HB stars: AV™, found for the GGCs, if a is as steep 
as <~ 0.3 mag/dex, globular clusters of different metallici- 
ties would be coeval; while if it is ~ 0.2, the most metal- 
rich globular clusters like 47 Tuc would be <~ 2 — 3 Gyr 
younger than the most metal-poor ones (Sandage 1993b, 
and references therein). 

Furthermore, the entire Population II distance scale 
would be affected if the Mv(RR) — [Fe/H] relation was 
found to be not universal and depending instead on the 
star formation history and the environment conditions of 
the host galaxy. 

It is thus important to derive the luminosity- 
metallicity relationship for RR Lyrae stars in different 
extragalactic systems, the LMC in particular. Walker 
(1992a) presented both photometry and metallicity esti- 
mates (based on the color magnitude diagrams and the 
properties of the RR Lyrae variables) for the LMC globu- 
lar clusters known to host RR Lyracs, which however are 
only a few (7) and cover a small range in metal abun- 
dance (-2.3 < [Fc/H] < -1.7, Walker 1992a, and refer- 
ences therein). Furthermore, these clusters are spread over 
a wide region of the sky, and some of them may well be 
at distances somewhat different from the average of the 
LMC: given the small numbers involved, uncertainties re- 
lated to the depth of the LMC may be as large as about 
0.1 mag. More recently, several thousands field RR Lyrae 
variables have been identified in the Clouds by surveys de- 
voted to the detection of microlensing events (MACHO: 
Alcock et al. 1996; OGLE: Udalski et al. 1997). However, 
metallicity estimates are available so far only for a handful 
of them (Alcock et al. 1996). If we combine the depth un- 
certainties with possible evolution off the ZAHB, at least 
100 variables are needed to obtain an accuracy of 0.05 
mag/dex in the slope of the luminosity-metallicity rela- 
tion. 

In recent years we began a study of the RR Lyrae 
variables in two fields close to the bar of the LMC. Our 
initial aim was to derive the mass-metallicity relation for 
the double mode pulsators (RRd), described in an ear- 
lier paper (Bragaglia et al. 2001). Accurate photometry 
for these fields in the Johnson BV and Cousins / bands 
was obtained using the Danish 1.5 m telescope, and spec- 



troscopy of the RRd variables was obtained using the 
EFOSC spectrograph at the 3.6 m ESO telescope. The 
analysis of the photometric data is described in sepa- 
rate papers (Clementini et al. 2003a, Di Fabrizio et al. 
2004). Adding spectroscopic observations, these data are 
well suited to a more extensive study of the properties of 
RR Lyrae stars in the LMC, in particular for determining 
their luminosity-metallicity relation, filling the gap left by 
previous studies. 

We were able to obtain adequate spectroscopic data 
for about 100 LMC RR Lyrae stars using FORS1 at VLT, 
and derived abundances using a variance of the Preston 
A S method (Preston 1959). An astrophysical discussion 
of the My (RR) — [Fe/H] relationship we derived, and of 
its impact on the distance to the LMC, has been already 
presented (Clementini et al. 2003a). Here, we give details 
on the acquisition and reduction of the spectroscopic data 
(Section 2), on the extraction of line indices (Section 3), 
and the derivation of the metallicities (Section 4). We 
then list the metal abundance of the individual variables, 
estimate the associated errors, discuss the properties of 
the sample (Section 5), compare the metal abundances of 
the a6-type RR Lyrae with those derived from analysis 
of the Fourier terms of the light curves (Section 6), and 
describe the derivation of the luminosity-metallicity rela- 
tion (Section 7). Finally, in Section 8 we comment on the 
metallicities of the short period Cepheids present in our 
sample. 

2. Observations and data reduction 

Spectroscopic data for the LMC targets and for HB 
stars in a number of calibrating globular clusters were 
collected with FORS1 (FOcal Reducer/low dispersion 
Spectrograph), mounted at UTl-Antu of the ESO Very 
Large Telescope (Paranal, Chile), on UT 21, 22, and 23 
December 2001. The first two nights were rather good (al- 
most photometric, with seeing varying between 0.5" and 
1.3", but generally less than 1"), and the last one slightly 
worse (light clouds present, and seeing varying between 
0.8" and 1.7", but generally larger than 1.2"). 

We used FORS1 in multi object (MOS) mode: up to 19 
spectra could be obtained in one exposure, using 19 pairs 
of movable slitlets, about 20" long, on a field of view (FoV) 
of 6.8 arcmin 2 . We selected the blue grism GRIS_600B, 
covering the 3450 - 5900 A wavelength range, with a dis- 
persion of 50 A mm -1 (R ~ 800), with slits 1" large. In 
this mode, each pixel is about 1.2 A. Only part of the 
FoV (about two thirds) was usable for each pointing, to 
cover the relevant wavelength range (i.e., ~ 3900 - 5100 
A) cointaining both the Can K and the hydrogen Balmer 
lines up to H/3. 

Pre-imaging frames were required, in conjunction with 
FIMS (the FORS Instrumental Mask Simulator package), 
to prepare all masks (42 in total: 33 on the LMC, 2 each 
on M68 and NGC1851, 1 on NGC3201, 4 on u Cen). In 
the case of the LMC, our original "Field A" and "Field 
B" (see Clementini et al. 2003a, and Di Fabrizio et al. 
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Fig. 1. F0RS1 LMC sub-field Al, Nord-East quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 

Fig. 2. FORS1 LMC sub-field A4, South-East quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 

Fig. 3. FORS1 LMC sub-field A2, Nord-West quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 

Fig. 4. FORS1 LMC sub-field A3, South- West quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 

2004), having a side of ~ 13 arcmin, were divided in 4 
sub-fields each to fit into the FORS usable FoV. We were 
able to fit into the 33 LMC masks about 100 of the 125 
RR Lyrae stars, and 3 of the 4 short period Cepheids we 
had detected and studied. 

Figures la-d, 2a-d show the eight LMC FORS1 sub- 
fields with the variables identified according to Di Fabrizio 
et al. (2004) identifiers. Centre of field coordinates are 
provided in Tabled 

Preparation of the observations was a delicate task, 
since we wanted to maximize the number of RR Lyrae's 
observed at/near minimum light for each pointing, as this 
minimizes errors in abundance derivations. We prepared a 
special software to help maximizing the efficiency of these 
observations. The nights were divided in time slots, con- 
sidering an overhead of 15 minutes for each pointing, and 
an exposure of 30 minutes on the LMC, and of 1-3 min- 
utes on the GC's. We then used our own positions, pe- 
riods, and epochs for the LMC RR Lyrae stars to derive 
for each time slot the maximum number of targets observ- 
able at minimum light in a single pointing, and filled the 
remaining slits with other RR Lyrae variables taken at 
random phase, and with clump stars. For the RR Lyrae 
stars in the calibration clusters M68, NGC 1851 and w Ccn 
we used the published light curves (from Walker 1998 for 
NGC 1851; from Clement et al. 1993, and Brocato et al. 
1994 for M68; from Kaluzny et al. 1997 for oj Cen), and 
re-derived ephemerides, to decide which RR Lyrae star to 
observe, choosing preferentially those at minimum light. 
For the NGC 3201 variables we used updated ephemerides 

Fig. 5. FORS1 LMC sub-field Bl, Nord-East quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 

Fig. 6. FORS1 LMC sub-field B4, South-East quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 



Fig. 7. FORS1 LMC sub-field B2, North- West quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 

Fig. 8. FORS1 LMC sub-field B3, South- West quadrant. 
Variables are marked by open circles and identified ac- 
cording to Di Fabrizio et al. (2004). 

Table 1. Centre of field coordinates of the eight LMC 



FORS1 sub-fields 


Field 


Q2000 


<52000 


Al 


5:23:36.47 


-70:34:23 


A2 


5:22:18.47 


-70:34:23 


A3 


5:22:18.47 


-70:40:43 


A4 


5:23:36.47 


-70:40:43 


Bl 


5:18:17.97 


-71:55:91 


132 


5:17:19.76 


-71:55:91 


B3 


5:17:19.76 


-71:66:51 


B4 


5:18:17.97 


-71:66:51 



and finding charts kindly provided to us in advance of pub- 
lication by A. Layden and A. Piersimoni, respectively. The 
efficiency of our selection is reflected in the large number of 
RR Lyrae that could be observed, since a simple division 
of the variable density for the available field and accept- 
able phase would have given about 50% less observations 
than we were actually able to perform. We obtained ad- 
equate spectroscopic data for 101 variables in the LMC, 
and observed also 9 RR Lyrae stars, 4 red and 1 blue 
HB stars in NGC 1851, 8 RR Lyrae's and 4 red HB's in 
NGC 3201, 13 RR Lyrae's, 7 blue and 1 red HB stars in 
M68, and 17 RR Lyrae variables, 1 Anomalous Cepheid 
(AC) and 1 Population II Cepheid (P2C) in oj Cen (see 
Tables 3, 4, 5 and 6). Multiple observations were obtained 
for about 50% of the LMC variables and for the RR Lyrae 
stars in ui Cen, in particular. Spectra were also obtained 
for 355 clump stars, their analysis is in progress. 

Exposure times were of 60 s for w Cen, 120 s for 
NGC3201, 180 s for M68 and NGC1851, and 1800 s (gener- 
ally, with a few 1600 and 2100 s exposures) for the LMC. 
For the latter, we had to compromise between signal to 
noise ratio (S/N) and time resolution, in order to ob- 
tain the highest S/N without phase smearing. The S/N 
varies, depending on exposure time, sky conditions, air- 
mass (about 1.4 to 1.7 in all nights), and on actual cen- 
tering of the star in the slit, but it is generally of about 
35 at 4700 A, and about 20 at 3950 A, as estimated by 
the pixel-to-pixel scatter. These values agree well with ex- 
pectations based on the Exposure Time Calculator when 
typical observing conditions are considered, thus showing 
that centering errors did not introduce significant light 
losses for most of the stars. 

Calibration frames (bias and flat field images, and 
wavelength calibration lamps) were obtained during day- 
time. 
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The MOS frames were reduced using standard routines 
in IRAF. 1 They were trimmed, corrected for bias and for 
the normalized flat field; special care was devoted to the 
flat-fielding procedure, since the CCD was read by four 
amplifiers, and a "normalization" was needed to eliminate 
jumps at the junctions. Up to 19 spectra were present in 
each frame, and were extracted with the optimal extrac- 
tion and automated cleaning options switched on. The sky 
contribution was subtracted making use of the slit length. 
Given the reasonable seeing conditions, contamination of 
targets from nearby stars was reduced to a minimum, ex- 
cept for a few objects. For each science mask, a HeCdHg 
lamp was acquired, and used to calibrate in wavelength 
the spectra, each one covering a different spectral range, 
depending on the target position, as usual in MOS obser- 
vations, but all comprising the Ca II K to H/3 lines. About 
10 lines of the calibration lamp were visible for each aper- 
ture, and the resulting dispersion solutions have r.m.s. of 
about 0.05 A. Further cleaning of cosmic rays hits and 
bad sky subtractions was done using the clipping option 
in the splot task; unfortunately a few spectra were unus- 
able for our procedure because of direct hits on the Ca II 
K line. Figure [5] shows examples of the obtained spec- 
tra. Details of the reduction procedures can be found in 
Taribello (2003). 

3. Line indices 

Metal abundances for the RR Lyrae variables can be ob- 
tained by comparing the strength of the Ca II K line with 
that of the H lines (Preston 1959). Preston used spectral 
types, best suited for photographic spectra; since we have 
CCD spectra, we prefer to use line indices. 

The spectra were measured using the program ROSA 
by Gratton (1989). The following steps were done on each 
spectrum. 

1. Geocentric radial velocities were measured by cross 
correlating the spectra with a suitable template (the 
RR Lyr variable 2 in M68, observed at high S/N). The 
zero point of the radial velocities was obtained by com- 
bining all data for M68, and comparing the resulting 
average velocity with that listed in Harris (1996). The 
error in the radial velocity measured from each indi- 
vidual spectrum, derived from the square average of 
the r.m.s. values for the 48 stars in our sample with 
multiple observations, is of 40 km s . This error in- 
cludes both the measure uncertainty, which is mainly 
due to centering errors of the stars on the slits, and 
the contribution due to variation of the radial velocity 
with phase during the pulsation. The error distribution 
of the 48 stars with multiple observations is gaussian, 
with only a few (6) outliers with errors larger than 70 

1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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MA) 

Fig. 9. Examples of spectra of stars in the bar of the LMC 
obtained using FORS1. The stars are identified according 
to Di Fabrizio et al. (2004). The seven upper plots are 
RR Lyrae variables observed at different phases. The two 
bottom plots are a short period Cepheid (star # 9604) and 
a much cooler clump star (star # 07376) , respectively. The 
spectra have been offset for clarity, and the main spectral 
lines are marked. 
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Wavelength (A) 

Fig. 10. Example of a spectrum of an RR Lyrae vari- 
able. Superposed are the spectral regions used to define 
the line indices. The solid lines define the spectral line 
bands; dashed lines define the comparison pseudocontin- 
uum bands. 
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Table 2. Definition of the spectral bands 




0.2 0.4 



Fig. 11. Correlation between the various indices of the 
H lines for the stars in M68 (filled squares), NGC 1851 
(open squares), and NGC3201 (crosses). Lines represents 
equality between indices. 



km s . The average velocity of these 48 stars, with- 
out applying phase corrections, is 261 ± 40 km s -1 (rms 
scatter), where 25 km s _1 is the internal error (since 
we have on average 2.4 observations per star) and 30 
km s _1 is the velocity dispersion. These radial veloci- 
ties are suited to study the kinematics of the old stellar 
component in the LMC, which however is beyond the 
purposes of the present work and will be discussed in 
a separate paper. 

2. All spectra were shifted to rest wavelength using the 
above measured geocentric radial velocities, and re- 
binned at costant wavelength step. 

3. In order to estimate line indices, instrumental fluxes 
in a few bands were measured on the spectra, by sim- 
ply integrating the spectra within given limits. These 
bands were centered on the Ca n K line, US, H/y, 
and H/3. For each line, we also defined two reference 



Index 


Min. Wavelength 


Max. Wavelength 


Weight 




(A) 


(A) 




Cn 


3910.0 


3925.7 


0.40 


K 


3925.7 


3941.7 




C12 


3941.7 


3953.0 


0.60 


C21 


4020.0 


4060.0 


0.50 


H,s 


4086.7 


4116.7 




C22 


4144.0 


4184.0 


0.50 


C31 


4260.5 


4300.5 


0.50 


1L, 


4325.5 


4355.5 




C32 


4380.5 


4420.5 


0.50 


C41 


4781.3 


4821.3 


0.29 


H/3 


4846.3 


4876.3 




C42 


4881.3 


4891.3 


0.71 



"pseudo-continuum" bands located roughly symmetri- 
cally on each side of the lines. The list of bands we 
used is given in Table [21 they are shown in Figure ITU1 

4. For each line i we constructed a line index 7^, defined 
as: 

h = 1 - Fi/[w a F(c a ) + w a F(c l2 )}, (1) 

where Fi is the instrumental flux measured for the 
band containing the line, F(cn) and F{ci2) are the 
fluxes measured in the reference continuum bands, and 
wn and Wi2 are weights that take into account the sep- 
aration in wavelength between the line and the com- 
parison bands. Note that wn + Wi2 = 1. By defining 
the line indices in this way we minimize the errors in 
the flux measurements due to the atmospheric disper- 
sion and centering on the slit, that cause part of the 
light from the star to fall out of the slit. 

5. The line indices for the three H lines are very well 
correlated to each other, except for a small trend of the 
Hp index, which is slightly larger for the cooler stars 
(see Figure ITT|) . We then defined an index < H > as 
simply the average of the three line indices measured 
for the hydrogen lines H^, H 7 , and H^. 

6. Finally, we defined K the line index for the Ca II K 
line. 

We did not apply any correction for the Ca II K inter- 
stellar absorption line, because it is expected to be small 
for the LMC and the calibrating cluster stars, that all 
have reddenings E(B — V) & 0.11 mag (see discussion in 
Gratton et al. 1986). 

Expected errors Oi in the indices may be estimated 
by differentating eq. (1). If we further assume that the 
average fluxes within each band are close to each other, 
and note that the sum of the weights of the comparison 
bands is equal to 1, we may write: 

a t » VA[l/^AX l +w tl /V^'i+w t 2/V^ r 2}/(S/N),(2) 

where A is the wavelength step, AA^, AAa, and AAj2 
are the widths of the bands containing the line i and the 
continuum comparison bands, and wn and wn are the 
weights attributed to them. 
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Table 3. Metallicity indices and metal abundances for the NGC 1851 stars 



Star 


Type 


HJD 




K 


Ha 


Y 


H s 


< H > 


M.I. 


[Re /HI 

L / J 


1 


ab 


2452266. 


.5229 





.325 


0.135 


0.130 


0.166 


0.144 


0.960 


-1.29 


5 


ab 


2452265. 


5278 


0. 


.427 


0.068 


0.069 


0.135 


0.091 


0.976 


-1.28 






2452266. 


5229 





.364 


0.085 


0.104 


0.122 


0.103 


0.803 


-1.42 


6 


ab 


2452266. 


5229 





.315 


0.157 


0.175 


0.178 


0.170 


1.162 


-1.13 


7 


ab 


2452265. 


.5278 


0. 


.116 


0.327 


0.350 


0.331 


0.336 


0.993 


-1.27 






2452266. 


.5229 





.326 


0.121 


0.128 


0.156 


0.135 


0.881 


-1.36 


11 


ab 


2452265. 


.5278 


0. 


.179 


0.259 


0.247 


0.247 


0.251 


0.964 


-1.29 






2452266. 


.5229 





.369 


0.134 


0.107 


0.143 


0.128 


1.043 


-1.23 


17 


ab 


2452266. 


.5229 


0. 


.226 


0.213 


0.218 


0.216 


0.216 


1.011 


-1.25 


21 


c 


2452266. 


.5229 





.182 


0.281 


0.264 


0.271 


0.272 


1.189 


-1.11 


22 


ab 


2452265. 


.5278 





.325 


0.133 


0.150 


0.193 


0.158 


1.101 


-1.18 


23 


c 


2452265. 


.5278 





.181 


0.287 


0.262 


0.270 


0.273 


1.186 


-1.11 






2452266. 


5229 





.126 


0.310 


0.293 


0.285 


0.296 


0.830 


-1.40 


99999 


bHB 


2452265. 


5278 





.150 


0.296 


0.279 


0.266 


0.280 


0.950 


-1.30 


129 


rHB 


2452266. 


5229 





.518 


0.043 


0.030 


0.125 


0.066 


1.132 


-1.15 


1031 


rHB 


2452266. 


.5229 





.407 


0.078 


0.072 


0.117 


0.089 


0.875 


-1.36 


1247 


rHB 


2452266. 


.5229 





.506 


0.016 


0.033 


0.102 


0.050 


0.950 


-1.30 


1278 


rHB 


2452266. 


.5229 





.534 


0.046 


0.036 


0.094 


0.058 


1.125 


-1.16 



Note: Star identifiers are from Sawyer- Hogg (1973) for the RR Lyrae stars, and from Walker (1992b) for the red HB stars 

(rHB). 



Table 4. Metallicity indices and metal abundances for the M68 stars 



Star 


Type 


HJD 




K 


Hp 


II- 


H s 


<H> 


M.I. 


[Fe/H] 


2 


ab 


2452264.8400 


0. 


.127 


0.158 


0.160 


0.167 


0. 


161 


-0.056 


-2.10 


8 


c 


2452265.8372 





081 


0.215 


0.167 


0.203 


0. 


.195 


-0.197 


-2.22 


9 


ab 


2452264.8400 


0. 


.075 


0.251 


0.228 


0.218 


0. 


232 


-0.060 


-2.11 






2452265.8372 


0. 


.124 


0.160 


0.158 


0.181 


0. 


.166 


-0.053 


-2.10 


12 


ab 


2452265.8372 


0. 


.107 


0.224 


0.201 


0.175 


0. 


.200 


0.010 


-2.05 


13 


c 


2452264.8400 


0. 


.119 


0.189 


0.168 


0.187 


0. 


.181 


-0.005 


-2.06 


17 


ab 


2452264.8400 


0. 


.112 


0.227 


0.219 


0.221 


0. 


.222 


0.184 


-1.91 


18 


c 


2452265.8372 


0. 


109 


0.219 


0.186 


0.210 


0. 


.205 


0.054 


-2.02 


20 


c 


2452265.8372 


0. 


102 


0.209 


0.207 


0.216 


0. 


.211 


0.038 


-2.03 


22 


ab 


2452265.8372 


0. 


.150 


0.130 


0.140 


0.176 


0. 


.149 


0.013 


-2.05 


24 


c 


2452264.8400 


0. 


129 


0.184 


0.169 


0.180 


0. 


.178 


0.044 


-2.03 


25 


ab 


2452264.8400 


0. 


.157 


0.156 


0.115 


0.134 


0. 


.135 


-0.013 


-2.07 


30 


ab 


2452264.8400 


0. 


164 


0.119 


0.119 


0.147 


0. 


128 


-0.011 


-2.07 


33 


c 


2452264.8400 


0. 


104 


0.213 


0.196 


0.206 


0. 


.205 


0.021 


-2.04 


86 


bHB 


2452264.8400 


0. 


014 


0.373 


0.372 


0.385 


0. 


.377 


0.030 


-2.04 


112 


bHB 


2452264.8400 


0. 


016 


0.392 


0.359 


0.359 


0. 


370 


0.026 


-2.04 






2452265.8372 


0. 


016 


0.364 


0.382 


0.355 


0. 


367 


0.016 


-2.05 


170 


bHB 


2452265.8372 


0. 


032 


0.320 


0.295 


0.297 


0. 


304 


-0.089 


-2.13 


340 


bHB 


2452265.8372 


0. 


.037 


0.331 


0.316 


0.292 


0. 


313 


0.012 


-2.05 


397 


bHB 


2452264.8400 


0. 


.006 


0.416 


0.399 


0.359 


0. 


391 


-0.023 


-2.08 


457 


rHB 


2452265.8372 


0. 


.176 


0.147 


0.105 


0.125 


0. 


126 


0.040 


-2.03 


464 


bHB 


2452265.8372 


0. 


043 


0.319 


0.309 


0.289 


0. 


305 


0.038 


-2.03 


535 


bHB 


2452264.8400 


0. 


Oil 


0.392 


0.379 


0.340 


0. 


.371 


-0.030 


-2.08 



Note: Star identifiers are from Sawyer-Hogg (1973) for the RR Lyrae stars, and from Walker (1994) for the blue and red HB 

stars (bHB and rHB, respectively). 



4. Metallicity calibration 

< H > is essentially a measure of the strength of the 
Hydrogen lines. Within the temperature range of interest, 

< H > is positively correlated with temperature, and it 
changes with phase since the temperature of a pulsating 
star changes with phase during the pulsation cycle. The 



run of < H > with phase for the globular cluster variables 
is shown in Figure 1121 In general, < H > is a function 
of the temperature T, gravity <?, and metal abundance 
[Fe/H], so that < H >=< H > (T,g, [Fe/H]). K is a mea- 
sure of the strength of the Ca II K line. In this range of 
temperatures, K decreases with temperature. Again, K is 
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Star 


Type 


HJD 




K 


Ha 


hi- 

i 


H s 


< H > 


M.I. 


[Re /HI 

L / J 


12 


ab 


2452266.8564 





.216 


0.171 


0.145 


0.179 


0.165 


0.506 


-1.65 


19 


ab 


2452266.8564 





.290 


0.128 


0.127 


0.169 


0.142 


0.745 


-1.46 


20 


ab 


2452266.8564 





.316 


0.106 


0.117 


0.168 


0.130 


0.790 


-1.43 


21 


ab 


2452266.8564 





.303 


0.116 


0.126 


0.169 


0.137 


0.782 


-1.43 


39 


ab 


2452266.8564 





.249 


0.133 


0.129 


0.176 


0.146 


0.552 


-1.62 


40 


ab 


2452266.8564 





.294 


0.136 


0.128 


0.144 


0.136 


0.719 


-1.48 


49 


ab 


2452266.8564 





.278 


0.158 


0.121 


0.150 


0.143 


0.689 


-1.51 


80 


ab 


2452266.8564 





.302 


0.121 


0.140 


0.156 


0.139 


0.787 


-1.43 


10 


rHB 


2452266.8564 





.423 


0.055 


0.067 


0.106 


0.076 


0.837 


-1.39 


263 


rHB 


2452266.8564 





.295 


0.122 


0.110 


0.150 


0.128 


0.660 


-1.53 


498 


rHB 


2452266.8564 





.461 


0.027 


0.043 


0.095 


0.055 


0.818 


-1.41 


2729 


rHB 


2452266.8564 





.345 


0.080 


0.079 


0.118 


0.092 


0.630 


-1.56 



Note: Star identifiers are from Sawyer-Hogg (1973) for the RR Lyrae stars, and from Rosenberg et al. (2000, available at 

http://dipastro.pd.astro.it/globulars/) for the red HB stars (rHB). 



Table 6. Metallicity indices and metal abundances for the variables in to Cen 



Star 


Typo 


HJD 




K 


Hp 


H-, 


H s 


< H > 


M.I. 


[Fe/H] 
This paper 


< [Fe/H\ > 
This paper 


fp /TT 

[Fe/H] 
BDE 


[Fe/H] 
caby 


<[Fe/H]> 
Others 


1 1 


c 


2452264 


.8479 





.121 


0.334 


0.328 


0.299 


0. 


321 


0.961 


- 1 .29 


-1 


.34 














2452264 


.8518 





.130 


0.324 


0.325 


0.306 


0. 


318 


1.041 


-1.23 


















2452265 


.8592 





.104 


0.341 


0.316 


0.281 


0. 


312 


0.720 


-1.48 














72 


P2C 


2452264 


.8479 





.354 


0.135 


0.129 


0.172 


0. 


145 


1.137 


-1.15 


-1 


.27 


-1.33 




-1 


.33 ± 0.30 






2452264 


.8518 





.360 


0.131 


0.127 


0.172 


0. 


143 


1.147 


-1.14 


















2452265 


.8592 





.298 


0.160 


0.142 


0.171 


0. 


158 


0.933 


-1.31 


















2452266 


.8462 





.345 


0.104 


0.097 


0.124 


0. 


108 


0.751 


-1.46 














73 


AC 


2452264 


.8479 





.187 


0.231 


0.211 


0.222 


0. 


221 


0.756 


-1.45 


-1 


.50 


-1.79 


-1.60 


-1 


.66 ± 0.17 






2452264 


.8518 





.191 


0.227 


0.205 


0.208 


0. 


213 


0.718 


-1.49 


















2452265 


.8592 





.171 


0.253 


0.242 


0.234 


0. 


243 


0.817 


-1.41 


















2452266 


.8462 





.139 


0.247 


0.230 


0.238 


0. 


238 


0.507 


-1.65 














74 


ab 


2452266 


.8462 





.265 


0.144 


0.127 


0.141 


0. 


137 


0.576 


-1.60 


-1 


.60 




-1.66 


-1 


.66 ± 0.20 


98 


c 


2452264 


.8479 





.130 


0.285 


0.261 


0.249 


0. 


265 


0.636 


-1.55 


-1 


.58 


-1.65 


-1.66 


-1 


.66 ± 0.17 






2452264 


.8518 





.127 


0.278 


0.256 


0.242 


0. 


259 


0.559 


-1.61 


















2452265 


.8592 





.097 


0.326 


0.296 


0.289 


0. 


303 


0.591 


-1.59 














99 


ab 


2452264 


.8479 





.225 


0.138 


0.118 


0.157 


0. 


138 


0.362 


-1.77 


-1 


.70 


-1.29 


-1.74 


-1 


.60 ± 0.17 






2452264 


.8518 





.238 


0.137 


0.111 


0.159 


0. 


136 


0.423 


-1.72 


















2452265 


.8592 





.262 


0.126 


0.136 


0.159 


0. 


140 


0.583 


-1.59 














100 


ab 


2452264 


.8479 





.153 


0.261 


0.257 


0.262 


0. 


260 


0.805 


-1.42 


-1 


.50 




-1.37 


-1 


.37 ± 0.20 






2452264 


.8518 





.148 


0.262 


0.258 


0.259 


0. 


260 


0.758 


-1.45 


















2452265 


.8592 





.145 


0.266 


0.266 


0.253 


0. 


262 


0.751 


-1.46 


















2452266 


.8462 





.107 


0.297 


0.271 


0.252 


0. 


273 


0.481 


-1.68 














124 


ab 


2452266 


.8462 





.267 


0.127 


0.098 


0.141 


0. 


122 


0.476 


-1.68 


-1 


.68 




-1.64 


-1 


.64 ± 0.20 


1(56 


ab 


2452265 


.8490 





.234 


0.133 


0.132 


0.158 


0. 


141 


0.435 


-1.71 


-1 


.71 


-0.69 


-1.67 


-1 


.37 ± 0.17 


168 


c 


2452265 


.8490 





.161 


0.259 


0.217 


0.225 


0. 


234 


0.653 


-1.54 


-1 


.54 


-1.45 


-1.33 


-1 


.37 ± 0.17 


160 


c 


2452265 


.8490 





.171 


0.271 


0.266 


0.261 


0. 


266 


1.023 


-1.24 


-1 


.24 




-1.64 


-1 


.64 ± 0.20 


170 


ab 


2452265 


.8490 





.258 


0.123 


0.132 


0.133 


0. 


129 


0.480 


-1.68 


-1 


.68 


-1.69 


-1.75 


-1 


.73 ± 0.17 


178 


c 


2452266 


.8368 





.163 


0.242 


0.183 


0.184 


0. 


203 


0.434 


-1.71 


-1 


.71 


-1.36 


-1.84 


-1 


.69 ± 0.17 


170 


ab 


2452266 


.8368 





.197 


0.139 


0.136 


0.164 


0. 


146 


0.265 


-1.85 


-1 


.85 


-1.14 


-1.78 


-1 


.58 ± 0.17 


182 


ab 


2452266 


.8368 





.208 


0.117 


0.125 


0.157 


0. 


133 


0.246 


-1.86 


-1 


.87 


-2.18 


-2.09 


-2 


.12 ± 0.17 


183 


ab 


2452266 


.8368 





.317 


0.162 


0.163 


0.193 


0. 


173 


1.200 


-1.10 


-1 


.10 


-0.75 


-1.49 


-1 


.26 ± 0.17 


184 


ab 


2452266 


.8368 





.333 


0.127 


0.111 


0.139 


0. 


126 


0.837 


-1.39 


-1 


.39 


-1.06 


-1.54 


-1 


.39 ± 0.17 


186 


c 


2452264 


.8479 





111 


0.347 


0.327 


0.313 


0. 


329 


0.905 


-1.34 


-1 


.22 














2452264 


.8518 





.117 


0.345 


0.326 


0.312 


0. 


327 


0.960 


-1.29 


















2452265 


.8592 





.144 


0.355 


0.318 


0.331 


0. 


335 


1.305 


-1.02 














187 


c 


2452266 


.8368 





.152 


0.246 


0.220 


0.223 


0. 


230 


0.556 


-1.62 


-1 


.62 




-1.59 


-1 


.59 ± 0.20 



Note: Star identifiers are from Kaluzny et al. (1997); the type classification for the Population II Ccphcid (P2C) and for the Anomalous Ccphcid 



(AC) is from Nemcc, Nemec, & Lutz 1994; BDE=Butlcr et al. (1978); caby= Rey et al. (2000); Others= average of BDE and Rey et al. (1978). 



a function of temperature, gravity, and metal abundance: 
K = K(T,g, [Fe/H]). If we now limit ourselves to stars on 
the HB 2 , we find that for these stars the gravity g is es- 
sentally a function of temperature and metal abundance. 
We may then simplify the relations for < H > and K and 

2 The three LMC short period Cepheids in our sample were 
all treated as being HB stars, we will come back to this point 
in Section 8. 



write < H >=< H > (T, [Fe/H]) and K = K(T, [Fe/H]). 
We further notice that the dependence of < H > on metal- 
licity is weak. We then expect that stars in a cluster (all 
having the same [Fe/H]) will define a nearly uniparamet- 
ric sequences in the < H > — K plane, with the sequences 
of different clusters being shifted according to metallicity. 
Mctallicities (for horizontal branch stars) can then be de- 
rived from the location of the stars in this plane, provided 
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ab-type 
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— I 1 1 1- 



c-type 




Phase 



Fig. 12. Run of the < H > index with phase for stars 
in the calibrating globular clusters. The upper panel is 
for fundamental mode pulsators {ab— type), and the lower 
panel for RR Lyrae stars pulsating in the first overtone 
(c-type). 



that suitable calibration sequences are available. This is 
the procedure adopted throughout this paper, following 
the recipe described below. 

First, we plotted the K index measured on each 
spectrum of the calibrating clusters M68 and NGC 1851 
against the < H > index (see Figure EI] and Table 3 and 
4) 3 . Note that Figure IT3l contains both variable and con- 
stant HB stars. We found that the RR Lyrae stars have 

< H > between 0.08 and 0.33: stars out of this range 
are non variables (those with < H > less than 0.08 are 
red horizontal branch stars, mainly observed in NGC1851; 
while those with < H > larger than 0.33 are blue hori- 
zontal branch stars, mainly observed in M68). RR Lyrae 
may be both fundamental mode (afr-type) and first over- 
tone (c-type) pulsators; when observed near minimum, ab- 
type RR Lyrae's have < H >^ 0.13, and c-type ones 

< H >^ 0.19. However, apart from segregation in the 
values of < H > when observed at minimum light, we 
found no systematic offset between the relations defined 
by variable and non variable stars whithin each cluster. 
In practice, we found that for each cluster there is a well 
defined relation between < H > and K, with a small scat- 
ter around it, although only for M68 and NGC 1851 the 
spread in < H > is large enough to define adequately the 
relation over a wide range of values of < H > . In fact for 
NGC 3201 we do not have blue HB stars and the variables 
were observed almost exclusively around minimum light, 



Fig. 13. Relation between the < H > and K spectral in- 
dices for stars in the calibration globular clusters: M68 
(squares), NGC 1851 (triangles). Different symbols are 
used for ab— type RR Lyraes (filled symbols), c— type RR 
Lyraes (open symbols) , and non variable HB stars (crosses 
with four and three arms, respectively). Superposed are 
the mean lines for M68 and NGC1851. A few of the vari- 
ables shown in the figure have multiple observations (see 
Tables 3 and 4). 

while, due to the spread in metallicity the w Cen vari- 
ables the relation between < H > and K for this cluster 
is very scattered. On the other hand, since the variables in 
M68 and NGC 1851 were observed when the stars where 
at different phases and may have different effective gravi- 
ties at the same temperature, we expected some thickness 
in these relations, since each star is expected to describe a 
loop in this diagram during its pulsation cycle. However, 
Figure El suggests that this thickness is small enough that 
can be neglected in our analysis. 

The mean relations drawn from the data in Figure HTH 

are: 



Ki = 0.3093 - 1.2815 < H > +1.3045 H 2 



(3) 



(valid for < H > between 0.12 and 0.40) for M68, that has 
[Fe/H]=— 2.06 according to the June 1999 update of Harris 
(1996) catalogue on Galactic globular clusters (available 
at http://www.physics.mcmaster.ca/Globular.html), and: 

t2 



K 2 = 0.6432 - 2.6043 < H > +3.0820 H l 



(4) 



NGC 3201 was not used here because the spectra of stars 
in this cluster all have similar < H > indeces 



(valid for < H > between 0.04 and 0.34) for NGC 1851, 
whose metallicity is [Fe/H]=— 1.26, again following Harris 
(1996). 

Metallicities for all other stars (both in the LMC and in 
NGC 3201 and wCen clusters) were then derived by linear 
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Fig. 14. Relations between the < H > and K spectral 
indices for variables in NGC3201 (upper panel) and u> Cen 
(lower panel). Superposed are the mean lines for M68 and 
NGC1851. 

interpolation/extrapolation between these two relations, 
entering the measured values of < H > and K indices. 
We defined a metallicity index: 

M.I. = {K- K x )j{K 2 - K x ) (5) 

where K is the Ca II K line index of the star and K\ , K 2 
are derived entering into equations (3) and (4) the < H > 
index measured for the star. The M.I. values derived by 
this procedure for the HB stars of the 4 calibrating clusters 
are listed in Columns 9 and 10 of Tables 3,4,5 and 6, and 
in Column 11 of Table 7 for the LMC variables. [Fe/H] 
abundances were then derived from the relation: 

[Fe/H] = [Fe/H] M68 + ([Fe/H] NGC1851 -[Fe/H] M68 )M.I.(6) 

Of course, these metallicities are less reliable when 
extrapolations outside the metallicity range defined by 
NGC 1851 and M68 are required. However, as it will be 




-2.5 -2 -1.5 -1 -0.5 



[Fe/H] (others) 

Fig. 15. Comparison between metal abundances from the 
present analysis and the average of the values determined 
by Butler et al. (1978) using the AS method, and by Rey 
et al. (2000) using the caby photometry, for RR Lyrae's 
in the globular cluster uj Cen. The solid line represents 
equality of the two values. 

discussed in Section 5 these extrapolations are neces- 
sary only for a few of the LMC variables in our sam- 
ple (see Figure 16). Metallicity indices and corresponding 
metal abundances obtained for the stars in the calibra- 
tion clusters using this procedure are given in Column 10 
of Tables 3, 4, 5, and 6. Metal abundances derived for the 
LMC variables are provided in Column 12 of Table and 
discussed in Section 5. 

As a first test of the accuracy of this abundance cal- 
ibration, we considered the case of NGC3201, for which 
8 variables were observed. All but one of the spectra for 
these stars were taken close to minimum light, i.e. over a 
rather small range of values of < H > (see upper panel 
of Figure 1141 and Table 5). The average abundance we ob- 
tained is [Fe/H]=— 1.49 ± 0.02, that nearly coincides with 
the value of [Fe/H]=-1.48 listed by Harris (1996). The 
star-to-star scatter of 0.08 dex is indeed very small, in 
agreement with the high S/N of these spectra. 

A further test is given by the abundances we obtain 
from the spectra of variables in lo Cen, a cluster that is 
known to have a large star-to-star scatter in the [Fe/H] val- 
ues for individual stars. Butler et al. (1978) obtained metal 
abundances for about 50 variables in this cluster, with 
average [Fe/H] =—1.43 for the ab— type variables (r.m.s. 
scatter of 0.43 dex), and [Fe/H]=— 1.72 for the c— type 
variables (r.m.s. scatter of 0.38 dex). Additional data were 
obtained by Rey et al. (2000) using the caby photometry. 
We obtained a total number of 36 spectra of 17 RR Lyrae 
stars, 1 Population II Cepheid and 1 Anomalous Cephcid 
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<H> 



Fig. 16. Relations between the < H > and K spectral 
indices for variables in the LMC. Superposed are the mean 
lines for M68 and NGC1851. The typical error bar is shown 
at the bottom-left corner. 

in this cluster (see Table 6 and lower panel of Figure 1X3)1 . 
For several stars we obtained spectra at different phases; 
by comparing the abundances obtained from these spec- 
tra, we obtain a standard error of [Fe/H] determination 
from a single spectrum of 0.11 dex, similar to the scat- 
ter obtained for NGC 3201. These [Fe/H] values may be 
compared with the [Fe/H] values obtained by averaging 
estimates by Butler et al. (1978) and Rey et al. (2000): 
these average values were obtained attributing errors of 
0.3 dex to the Butler et al.'s [Fe/H] values, and of 0.2 dex 
to those of Rey et al., in agreement with the typical er- 
rors estimated in the original papers. This comparison is 
shown in Figure 1151 The agreement is good: the average 
difference is 0.01 ± 0.04 dex, with a star-to-star scatter of 
the residuals of 0.18 dex, in good agreement with the error 
bars both in our and other determinations. 

5. Metal abundances of the LMC variables 

Figure rj^| shows the run of the K spectral index as a func- 
tion of the < H > for the variables in the two fields of 
the LMC. Most of the stars lie between the ridge lines for 
M68 ([Fe/H]=-2.06) and NGC1851 QFe/H]=-1.26), with 
only a few stars above the latter line (i.e. more metal-rich 
than this cluster). We obtained metallicity determinations 
for 101 LMC variables from the analysis of 168 spectra. 
According to Di Fabrizio et al. (2004) 98 of these stars 
are RR Lyrae, and 3 are short period Cepheids. Table 
lists the metal abundances for the stars in the two fields 
of the LMC. Column 1 gives the variable identification 
with the first two characters indicating the FORS1 sub- 



field where the star is located (see Figures la-d, 2a-d), and 
the remaining digits giving Di Fabrizio et al. (2004) iden- 
tification number. Column 2 gives the variable type (ab: 
fundamental mode; c: first overtone; d: double-mode pul- 
sator; SPC: Short Period Cepheid). Columns 3 and 4 give 
the mean observed and dereddened V magnitudes taken 
from Di Fabrizio et al. (2004). According to Clementini 
et al. (2003a) reddenings of E(B - V)=0.116 and 0.086 
mag were adopted for the stars in field A and B, respec- 
tively. Column 5 gives the Heliocentric Julian Day (HJD) 
of observation at mid exposure. Columns 6 to 9 give the 
measured line indices, and Column 10 gives the value of 

< H > (the average for the three H line indices). Column 
11 gives the metallicity index (M.I.) and Column 12 gives 
the metallicity obtained for each spectrum from our cali- 
bration of the line indices. Tabulated errors are obtained 
as detailed below. Finally, Column 13 gives the final value 
of the metallicity for each variable obtained by a weighted 
average of the metallicity determinations from individual 
spectra. These metallicities are tied to Harris (1996) metal 
abundance for the Galactic globular clusters NGC 1851, 
NGC3201, and M68 ([Fe/H]=-1.26, -1.48, and -2.06, 
respectively, to compare with —1.33, —1.56, and —2.09 of 
Zinn & West 1984 metallicity scale), that we used as cal- 
ibrators. Thus our metallicities are on a metallicity scale 
that, on average, is 0.06 dex more metal rich than the Zinn 
& West one 4 . 

We note that all stars have [Fe/H]< —0.9, except A3- 
05148, that is flagged as a blend by Di Fabrizio et al. 
(this is the only star with this label in our sample), A3- 
02119 and Al-10360, c— type variables for which we ob- 
tain [Fe/H]=-0.79±0.26 and -0.27 ±0.24, respectively. 
It should also be noted that abundances for stars whose 
metallicity is larger than [Fe/H] =—1.26 (the metallicity 
of NGC1851) are actually obtained by extrapolation, and 
are then much more uncertain: this applies in particular 
to the stars with [Fe/H]> —1.0. On the whole we think 
that there is little evidence for RR Lyrae stars more metal 
rich than [Fe/H]= —1 in the LMC, from our data. 

To estimate observational errors in the [Fe/H] values, 
we have compared results obtained from multiple obser- 
vations of the same star. Since the ridge lines for the two 
comparison clusters tend to converge at large values of 

< H > (i.e., high temperatures, that is spectra taken far 
from minimum light), we expect that the accuracy of the 
metallicities is a function of the phase at which the spec- 
tra were taken. We have then divided the stars into three 
groups according to the average strength of the < H > pa- 
rameter. For each group we estimated the quadratic mean 
of the r.m.s., weighting values for individual stars accord- 
ing to the number of observations; results are shown in 



4 There is no clear offset between the Harris (1996) and Zinn 
& West (1984) metallicity scales, the mean offset being only 
0.02±0.01 dex (113 clusters), once a few clearly discrepant ob- 
jects are eliminated. We consider the 0.06 dex difference found 
for the three calibrating clusters of the present analysis as a 
possible zero point error in our abundance scale. 
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Table 7. Metallicities for LMC RR Lyrae stars and short period Cephcids 
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typo 


V 


V a 


HJD 




K 


H s 


H-, 


Hp 


< H > 


M.I. 




[Fe/Hj 




<[Fc/H]> 
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06332 


c 
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2452264. 
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0. 
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0. 
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1 
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— 
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-0 
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d 
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0. 
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0. 
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27 


-1 
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27 
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c 
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18 
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0. 
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0. 
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-1 
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■07247 


ab 
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0. 
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.848 
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-1 
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10 














2452265. 
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0. 
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.129 





.121 
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.806 


— 1 


.42 


± 


0. 


11 
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ab 


19.435 


19 


103 
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6233 
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0. 
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.271 


1 


.098 


— 1 


.18 


± 





26 


-1 


.28 


± 





00 
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0. 
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± 
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0. 
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ab 
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18 
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2452265. 
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0. 
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± 





28 
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28 
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ab 
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18 
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0. 
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± 


0. 


11 
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0. 
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0. 
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0. 
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0. 
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± 


0. 


12 
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.83 


± 


0. 


12 
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d 
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0. 
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0. 
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— 1 


.50 
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c 
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0. 
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0. 
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22 
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22 
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ab 
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0. 
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.157 
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.66 


± 


0. 


14 


-1 
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± 


0. 


14 
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ab 
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0. 
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.986 
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18 


-1 
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.184 





.198 


1. 
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± 
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17 
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± 
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± 
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10 
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ab 
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19.283 
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0. 
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0. 
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1. 
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Table |H1 As expected, the r.m.s. scatter increases with in- 
creasing strength of the H lines. We then made a linear fit 
throughout these data, and assumed that the error in the 
metallicity determination is given by: 

cr[Fe/H] = 1.027 < H > -0.017. (7) 

These empirical estimates of the errors in our metallic- 
ities can be compared with the values expected from the 
quality of our spectra (that we remind are within expec- 
tations on the basis of the Exposure Time Calculator). 
Errors are mainly due to the K index, the contribution 
by errors in the < H > index (evaluated to be 0.014 by 
comparing the values obtained from different lines) being 
smaller, although not entirely negligible. By applying eq. 
(2), the expected error in the K index is 0.029 for a typ- 
ical spectrum in the LMC. This value may be compared 
with that obtained from multiple spectra of the same star 
taken at very similar phase. Using eight independent pairs 
of observations for which the < H > index differs less than 
0.01, the error in individual K measurements is 0.033, in 



good agreement with expectations. If we now take into 
account the sensitivity of metallicity to variations of the 
K index (as a function of < H >), that can be obtained 
by combining eqs. (3) and (4), and the small contribution 
given by the errors in < H > (provided by the slopes of 
the eqs. (3) and (4)), we may estimate the predicted errors 
for typical RR Lyrae stars of the LMC: the expected val- 
ues for the three bins of Table 00 are 0.157, 0.189, and 
0.226 dex, respectively. These values, listed in the last 
column of Table agree well with the observed errors, 
confirming that the accuracy of our abundance determi- 
nations is within expectations. 

Multiple observations of the same variable were then 
combined, weighting individual determinations according 
to the error estimated in this way thus producing the av- 
erage values listed in Column 13 of Table [7| 

The metallicity distribution of the RR Lyrae stars in 
our two fields of the LMC is given in Table Eland shown 
in Figure El The average metallicity of the sample is 
<[Fe/H>= -1.48 ± 0.03 ± 0.06 dex, where the last error 
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Table 7. - continued - 
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bar accounts for possible zero point errors in our calibra- 
tion. The star-to-star scatter (r.m.s=0.29 dex) is clearly- 
larger than the mean quadratic error of 0.17 dex of the 
abundances for individual stars 5 . By quadratically sub- 
tracting the two values, the intrinsic star-to-star scatter 
is 0.23 dex. There is no clear systematic offset between 
field A and B: the average abundances are —1.45 ± 0.04 
and —1.52 ± 0.04 dex, respectively. The distribution of 
Figure HTI aDpears to be somewhat skewed, with more stars 
more metal-poor than the average value than stars more 
metal rich than this value. 

If we divide the variables according to the Bailey types, 
we find average metallicities of [Fe/H]=— 1.51 ± 0.03 dex 
for ab-type RR Lyrae's, [Fe/H]=-1.36 ± 0.07 dex for 
c— type ones, and finally [Fe/H]=— 1.57 ± 0.09 dex for 
double pulsators. Differences between values for different 
types are quite small; the somewhat larger value found 



5 This error is smaller than the values quoted in Table |S| 
because on average about 1.6 observations were obtained for 
each star 



Table 8. Errors as a function of < H > for the LMC stars 



< H > bin n. stars < H > r.m.s.([Fe/H]) Expected 



< 0.17 


16 


0.151 


0.138 


0.157 


0.17-0.22 


18 


0.203 


0.193 


0.189 


> 0.22 


14 


0.259 


0.249 


0.226 



for c— type variables may be due to a few stars with larger 
error bars, for which we found metallicities [Fe/H]> — 1. 
Note that since on average c— type variables have higher 
temperatures, they have weaker Ca n K lines, making the 
errors of the metal abundance determinations larger. 

6. Comparison with abundances from light curve 
Fourier analysis 

Metallicities from the parameters of the Fourier decom- 
position of the V light curves were derived for 29 RRab's 
in our photometric sample, 7 of which do not have spec- 
troscopic metal abundances, applying Jurcsik & Kovacs 
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Table 7. (- continued -) 



Star 


typo 


V 




H.JD 


K 


H 5 




Hp 


< H > 


M.I. 


TF — /ttI 

[Fe/HJ 




<lFc/H]> 


Bl-05952 


SPC 


18.459 


18.192 


2452265.7069 


0.378 


0.044 


0.050 


0.125 





.073 


0.619 


-1.56 ± 


.06 


-1.59 : ±0.03 










2452265.8201 


0.389 


0.060 


0.048 


0.096 


0. 


068 


0.642 


-1.55 ± 


.05 














2452266.5993 


0.419 


0.005 


-0.001 


0.088 


0. 


031 


0.502 


-1.66 ± 


.05 














2452266.6840 


0.209 


0.207 


0.170 


0.193 


0. 


190 


0.657 


-1.53 ± 


.18 






Bl-06164 


c 


19.057 


18.818 


2452265.8201 


0.111 


0.234 


0.217 


0.241 


0. 


.230 


0.225 


-1.88 ± 


.22 


— 1 


.79 ± 0.15 










2452266.5993 


0.157 


0.198 


0.207 


0.224 


0. 


210 


0.436 


-1.71 ± 


.20 






Bl-06957 


c 


19.197 


18.958 


2452265.7069 


0.223 


0.216 


0.167 


0.176 


0. 


187 


0.719 


-1.48 ± 


.18 


-1 


.48 ± 0.18 


Bl-07063 


ab 


19.629 


19.390 


2452265.8201 


0.273 


0.168 


0.131 


0.150 


{) 


.150 


0.719 


-1.49 ± 


.14 


— 1 


.49 ± 0.14 


Bl-07064 


c 


19.122 


18.883 


2452266.5993 


0.100 


0.207 


0.198 


0.234 
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0.032 


-2.03 ± 


.20 


— 1 


.96 ± 0.16 










2452266.6840 


0.092 


0.277 


0.251 


0.260 


0. 


263 


0.272 


-1.84 ± 


.25 
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19.043 


18.804 
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0.126 


0.240 


0.221 


0.230 





.230 


0.341 


-1.79 ± 


.22 


— 1 


.64 ± 0.10 










2452265.8201 


0.175 


0.233 


0.214 


0.226 


0. 


224 


0.696 


-1.50 ± 


.21 














2452266.5993 


0.165 


0.192 


0.147 


0.194 


0. 


178 


0.273 


-1.84 ± 


.17 














2452266.6840 


0.179 


0.257 


0.321 


0.236 


0. 


272 


1.146 


-1.14 ± 


.26 
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ab 


19.078 


18.839 


2452265.7069 


0.164 


0.149 


0.126 


0.129 


0. 


135 


0.017 


-2.05 ± 


12 


—2 


.05 ± 0.12 
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ab 


19.426 


19.187 


2452265.8201 


0.278 


0.156 


0.168 


0.196 


0. 


173 


0.957 


-1.29 ± 


.16 


-1 


.34 ± 0.11 










2452266.5993 


0.285 


0.121 


0.177 


0.174 


0. 


157 
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-1.38 ± 


.14 
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ab 


19.385 


19.146 


2452265.7069 


0.250 


0.181 


0.119 


0.169 


0. 


156 


0.641 


-1.55 ± 


.14 


-1 


.43 ± 0.09 
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0.204 


0.181 


0.185 


0.206 


0. 


191 


0.625 


-1.56 ± 


.18 
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0.237 
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0.242 


0.194 
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217 


1.103 


-1.18 ± 


.21 
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19.126 
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0.242 


0.152 


0.165 


0.203 


0. 


.173 


0.730 


-1.48 ± 


.16 


— 1 


.31 ± 0.09 
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0.329 


0.183 


0.161 


0.162 


0. 
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-1.07 ± 
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0.139 


0.150 
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.20 ± 0.08 
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0.317 
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0.140 


0.156 


0.178 
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— 1.16 ± 


.15 
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0.243 


0.216 


0.229 
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0.771 


-1.44 ± 
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0.132 


0.112 


0.131 


0.158 





.133 
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.12 


—2 


,12 ± 0.11 










2452265.5716 


0.092 


0.275 


0.264 


0.261 


0. 


266 


0.290 


-1.83 ± 


.26 






B2-06020 


ab 


Shift 
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0.123 


0.145 


0. 
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-1 


.86 ± 0.11 
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0.220 


0.182 


0.151 
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0.057 
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0.209 
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-1 


.52 ± 0.16 
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ab 
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0.274 


0.272 


0.280 


0. 


275 


1.083 


-1.19 ± 


.27 


-1 


.19 ± 0.27 
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18.968 


2452266.7669 
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0.120 


0.038 
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0. 


.108 
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-1 


.48 ± 0.09 
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0.169 
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0.238 
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19.014 
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0.170 


0.278 


0.244 


0.211 


0. 


.244 


0.821 


-1.40 ± 


.23 


— 1 


.20 ± 0.13 










2452265.5716 


0.314 


0.159 


0.169 


0.191 
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173 


1.184 


-1.11 ± 
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2452265.5716 
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0.112 


0.104 


0.165 
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.127 


0.600 


-1.58 ± 
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— 1 


.58 ± 0.11 


B2-07648 
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2452265.5716 


0.127 
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0.291 


0.305 





.308 
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-1.32 ± 


.30 
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.32 ± 0.30 
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0.070 
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.26 ± 0.12 
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2452266.8216 
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0.165 
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.122 


0.450 


-1.70 ± 
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— 1 


.70 ± 0.11 
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0.140 
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0.268 
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0.217 
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.90 ± 0.09 
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0.199 


0.121 


0.133 


0.143 
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-1.91 ± 
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.189 
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.65 ± 0.13 
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0.209 
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.254 
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0.133 
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0.935 


-1.31 ± 


.29 














2452266.8216 
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0.285 
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19.042 


2452265.7350 


0.136 


0.303 
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0.268 





.286 
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.37 ± 0.28 
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19.173 


18.934 
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0.121 


0.284 


0.282 


0.274 
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0.668 


-1.53 ± 


.27 


— 1 


.Oo ± 0.2 / 
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19.514 


19.275 
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0.124 
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au 


19.286 


19.047 


z4t0ZZ04.0OO± 


U. 1U4 




n o^o 
u.zoy 


U.zDo 





273 


n 07Q 
U.z / o 


— 1. fU ± U 


Of. 


-1 


.52 ± 0.18 










2452264.5656 


0.163 


0.275 


0.241 


0.257 


0. 


258 


0.258 


-1.36 ± 


.25 






B4-02811 


c 


19.384 


19.145 


2452264.5351 


0.205 


0.245 


0.219 


0.243 


0. 


.236 


0.236 


-1.23 ± 


.23 


-1 


.24 ± 0.14 










2452264.5656 


0.254 


0.171 


0.185 


0.235 


0. 


197 


0.197 


-1.24 ± 


.19 






B4-03033 


ab 


Incomplete 


2452264.5351 


0.219 


0.232 


0.212 


0.215 





.220 


1.000 


-1.26 ± 


.21 


-1 


.36 ± 0.15 










2452264.5656 


0.168 


0.266 


0.231 


0.208 


0. 


235 


0.235 


-1.48 ± 


.22 






B4-04749 


c 


19.314 


19.075 


2452264.5351 


0.246 


0.175 


0.170 


0.176 





.174 


0.174 


-1.45 ± 


.16 


-1 


.40 ± 0.13 










2452264.5656 


0.220 


0.198 


0.223 


0.218 





213 


0.213 


-1.31 ± 


.20 






B4-10811 


ab 


19.431 


19.192 


2452264.5656 


0.201 


0.231 


0.187 


0.223 


0. 


214 


0.214 


-1.42 ± 


.20 


-1 


.42 ± 0.20 



Note: V and Vq are mean apparent and dereddened intensity averaged magnitudes taken from Di Fabrizio et al. (2004). 



Shift means that there is an offset between light curves obtained in 1999 and in 2001; Incomplete means that the light 
curve is incomplete, so that no meaningful average magnitudes could be obtained; Blazhko means that the light curve 
shows evidence of a significant Blazhko effect (Blazhko 1907); Blend means that the stellar image appears blended in 
some frames; The metal abundance of Bl-05952 is uncertain (see Section 8). 



(1996) and Kovacs & Walker (2001) techniques (see Di 
Fabrizio et al. 2004, for details). 

These metal abundances are on the metallicity scale 
defined by Jurcsik (1995) which is, on average, about 0.2 
dex more metal rich than Harris (1996) scale. The aver- 



age metallicity of this subsamplc of stars is: [Fe/H]=— 1.27 
(cr=0.35, 29 stars), in good agreement with the average 
values derived from the spectroscopic analyses, once dif- 
ferences between metallicity scales are properly taken into 
account. 
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Fig. 17. Metallicity distribution for RR Lyrae stars in the 
LMC. The bottom panel shows the distribution of all vari- 
ables; the other three panels show the distribution of ab~ , 
c— and d— type RR Lyrae stars, respectively. 

Table 9. Metallicity distribution and average magnitudes 
of the RR Lyrae stars in the LMC 



[Fe/H] 


bin 








n. stars 


< [Fc/H] > 


< Vo > 


[Fe/H] < - 


-2.1 






1 


-2.12 


18.88 


-2.1 < 


[Fe/H] 


< - 


2 





3 


-2.04 


18.91 ±0.05 


-2.0 < 


[Fe/H] 


< - 


1 


9 


3 


-1.93 


19.15 ±0.15 


-1.9 < 


[Fe/H] 


< - 


1 


8 


4 


-1.84 


19.12 ± 0.08 


-1.8 < 


[Fe/H] 


< - 


-1 


7 


7 


-1.73 


19.04 ± 0.08 


-1.7 < 


[Fe/H] 


< - 


-1 


6 


12 


-1.64 


19.02 ±0.05 


-1.6 < 


[Fe/H] 


< - 


1 


5 


15 


-1.54 


19.10 ±0.03 


-1.5 < 


[Fe/H] 


< - 


1 


4 


21 


-1.45 


19.09 ± 0.04 


-1.4 < 


[Fe/H] 


< - 


1 


3 


18 


-1.35 


19.07 ±0.03 


-1.3 < 


[Fe/H] 


< - 


-1 


2 


6 


-1.24 


19.20 ±0.07 


-1.2 < 


[Fe/H] 


< - 


-1 


1 


1 


-1.19 


19.01 


-1.1 < 


[Fe/H] 


< - 


-1 












-1.0 < 


[Fe/H] 


< - 





9 


4 


-0.92 


19.30 ±0.10 


-0.9 < 


[Fe/H] 


< - 





8 









-0.8 < 


[Fe/H] 


< - 





7 


1 


-0.79 


19.33 


-0.7 < 


[Fe/H] 








2 


-0.37 





The star-to-star comparison between photometric and 
spectroscopic abundances is shown in Figure 1181 On av- 
erage, the difference is 0.30 ± 0.07 dex, with photometric 
abundances being larger as expected. However, there are 
three stars (A10214, A26525, and B3-01408) that have 
large residuals between photometric and spectroscopic 
abundances. The light curve of one of these stars (A10214) 
does not fully satisfy the completness and regularity cri- 
teria for a reliable application of Jurcsik & Kovacs (1996) 
method (the star has a large D m value, see Di Fabrizio et 
al. 2004) and, for all these stars the reddening E{B — V) 




-1 .5 
[Fe/H] (phot) 

Fig. 18. Comparison between metal abundances for LMC 
RR Lyrae's obtained from spectroscopy and from analysis 
of the terms of a Fourier decomposition of the light curves 
following Jurcsik & Kovacs (1996) and Kovacs & Walker 
(2001) technique. Open symbols mark stars for which dis- 
crepant reddening values were derived using Sturch (1966) 
method. The solid line is the mean offset between the two 
_ metallicity sets. 

derived from the Fourier decomposition of the light curve 
is quite discrepant with respect to the average of the other 
stars. This suggests that the Fourier analysis for these 
stars may be unreliable. If they are dropped, the average 
difference between the two sets of metal abundances is 
0.28±0.05 dex, in reasonable agreement with the expected 
offset between the Jurcsik (1995) and Harris (1996) metal- 
licity scales. The r.m.s of the residuals (0.24 dex), agrees 
well with observational errors of both metallicity deter- 
minations (0.17 dex from spectroscopy; ~ 0.2 dex from 
photometry) . 

7. The luminosity - metallicity relation for RR 
Lyrae stars 

As discussed in the Introduction the metallicity [Fe/H] 
influences the absolute magnitude of the RR Lyrae stars 
My (RR) through a relation generally assumed to be of lin- 
ear form: My(RR)=ax [Fe/H]+/3. Rather large uncertain- 
ties exist on the slope a of this relation whose derivation 
was the ultimate goal of the present spectroscopic study. 

The LMC bar is the ideal place for defining the slope 
of the My (RR)- [Fe/H] relation since it contains a very nu- 
merous population of RR Lyrae stars, spanning more than 
1.5 dex in metallicity, and the variables are all at the same 
distance from us, hence we do not have to worry about the 
absolute magnitudes, i.e. about model assumptions. There 
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Table 10. Metallicity bins used to calculate the 



luminosity-metallicity relation 


of the LMC RR Lyrae stars 

v 


1B.0 




[Fc/H] bin 




n. stars < [Fc/H] > 


< V > 


1 8 5 


; 


[Fe/H] < - 


-1.8 


9 


-1.96 ±0.06 


18.99 ± 0.05 






-1.8 < [Fe/H] < - 


-1.6 


18 


-1.67 ±0.04 


19.01 ± 0.04 


19.0 




-1.6 < [Fe/H] < - 


-1.4 


32 


-1.49 ±0.03 


19.07 ±0.03 






-1.4 < [Fe/H] < - 


-1.2 


21 


-1.32 ±0.04 


19.09 ±0.03 


19.5 




-1.2 < [Fe/H] 




7 


-0.97 ±0.07 


19.21 ± 0.06 







V 0.214([Fe/II] + 1.5) + 19.064 



still remains to consider the intrinsic spread in the lumi- 
nosities of the RR Lyrae stars due to their evolution off 
the ZAHB. To minimize this effect a fairly large number 
of variables is needed: we estimated that about 100 stars 
are required to reduce the error bar in a down to better 
than 0.05 mag/dex. 

We have combined the individual metal abundances of 
the RR Lyrae stars derived from our spectroscopic study 
(see Column 13 of Table[7J) with the corresponding dered- 
dened mean apparent V magnitudes derived from our pho- 
tometric work (see Column 4 of Table |7J) to determine the 
slope of the luminosity-metallicity relation for RR Lyrae 
stars with accuracy of 0.05 mag/dex. For 11 of the RR 
Lyrae variables we have analyzed spectroscopically, the 
light curve is either incomplete, or with systematic shifts 
between the 1999 and the 2001 photometry, or the star im- 
age appears to be blended (see Columns 3-4 of Table [TJ), 
so that accurate values of the average apparent magni- 
tude could not be obtained for these stars. Our luminosity- 
metallicity relation is thus based on the remaining 87 stars 
which cover the metallicity range from [Fe/H]=— 0.79 to 
—2.12. We have divided the whole metallicity range into 
5 bins, to have a reasonable number of objects in each of 
them, and computed both the average metal abundance 
and the average dereddened apparent luminosity Vq and 
the corresponding rms errors for each bin, these values are 
given in Table ITUl 

From a least square fit of these average values weighted 
by the errors in both variables we derive the following rela- 
tion between the apparent luminosity and the metallicity 
of the RR Lyrae stars in our sample: 



< Vq >= 0.214(±0.047)([Fe/H]±1.5)±19.064(±0.017)(8) 

where the error in the slope was evaluated via Monte Carlo 
simulations. The relation is shown in Figure 19, where 
we distinguish between single and double-mode RR Lyrae 
stars, indicated by open and filled squares respectively. 
We note that the LMC RRd's are systematically offset 
to slightly higher luminosities, thus suggesting that these 
stars may be more evolved than their single-mode pulsator 
counterparts. 

A complete discussion of the impact of the 
new luminosity-metallicity relation has been given in 
Clementini et al. (2003a); here we only recall that the 
mild slope we derive agrees very well both with the re- 
sults by Rich et al. (2001) based on the HB luminosity of 
19 globular clusters in M31 (AM v (RR)/A[Fe/H] = 0.22 
mag/dex), and with the 0.20±0.04 slope found by the 



-2 -1.5 -1 

[Fe/H] 



Fig. 19. Run of the dereddened mean magnitude of the 
RR Lyrae stars in our two fields of the LMC as a function 
of the metallicity [Fe/H]. Filled symbols are double- mode 
pulsators; open symbols ab- and c-type RR Lyrae stars. 

Baade-Wesselink analysis of Milky Way field RR Lyrae 
stars (Fernley et al. 1998). We also note that Figure 19 
does not show clear evidence for a break in the slope 
around [Fe/H]=— 1.5 as proposed by the theoretical mod- 
els. This seems to suggest that the luminosity-metallicity 
relation for HB stars is universal, at least in M31, the 
LMC and the Milky Way. However, we caution the reader 
that there are only 4 stars with metallicity larger than 
[Fe/H]=— 1 in our sample. A much larger number of stars 
with metallicities higher than [Fe/H] = — 1.0 would be nec- 
essary to definitely asses whether the My (RR)- [Fe/H] re- 
lation actually breaks at [Fe/H]=— 1.5. 

8. The average metallicity of the short period 
Cepheids 

Our LMC sample includes 3 variables (star #9604 and 
#10320 in Field A, and star #5952 in field B) that have 
periods (0.29< P < 0.63 days) in the range typical of the 
RR Lyrae stars but are from 0.5 to about 0.9 mag brighter 
than the average luminosity of the HB in these fields (see 
left panels of Figures 4 and 5 of Clementini et al. 2003a). 
The classification of these variables posed some difficulty 
since, as discussed in Di Fabrizio et al. (2004), their av- 
erage luminosities, periods and amplitudes are consistent 
either with those of the Anomalous Cepheids (ACs) com- 
monly found in dwarf Spheroidal galaxies (e.g. Pritzl et 
al. 2002, and references therein), or of the low luminosity 
(LL) Cepheids (Clementini et al. 2003b) and of the short 
period Classical Cepheids found in a number of dwarf 
Irregular galaxies (Smith et al. 1992, Gallart et al. 1999, 
2003, Dolphin et al. 2002). 

Knowlewdge of the metallicy may help to classify these 
variables since there is a limiting metallicity above which 
no ACs should be generated (Bono et al. 1997, Marconi 
et al. 2004) 6 , and ACs are expected to have low metal 
abundances, similar to or lower than the metallicity of 
the oldest populations in the host system, while short pe- 



6 This limit should be [Fe/H]~ —1.7 for variables around 
- 1.3M and [Fe/H]~ -2.3 for variables around ~ 1.8M . 
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\ (A) 

Fig. 20. Spectra of the variable star #73, one of the two 
suspected Anomalous Cepheids in uj Cen. For each spec- 
trum we give the phase along the pulsation cycle; the mini- 
mum light corresponds to the phase interval 0.40< <fi m i n < 
0.52 (bottom spectrum). 
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Fig. 21. Spectra of the 3 LMC SPC in our sample. For 
each spectrum we give the corresponding phase along the 
pulsation cycle. The interval of minimum light is 0.77< 
4>min < 0.85 for star A2-9604, 0.50 < <f> mirl < 0.68 for star 
A2-10320, and 0.77 < <j) min < 0.83 for star Bl-5952. 

riod Classical Cepheids should have metallicities similar 
to those of the Population I component in the system. 

Two suspected ACs are known in u> Cen (Nemec et al. 
1994, Kaluzny et al. 1997), a cluster known to span a wide 
range in metallicity with at least three separate enrich- 
ment peaks (Norris, Freeman, & Mighell 1996, Suntzeff & 
Kraft 1996, Pancino et al. 2002) and suggested to possi- 
bly be the remnant of a disrupted dwarf galaxy. One of 
them is among the stars in our sample (variable #73, see 
Table 6). The spectra of this star are shown in Figure I2U1 
those of the 3 LMC short period Cepheids are shown in 
Figure ED 

In the derivation of the spectroscopic abundances both 
the Anomalous Cepheid in cu Cen (variable #73), and the 
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Fig. 22. Run of < H > with the dereddened B — V colour 
for the LMC variables observed around minimum light. 
Different symbols are used for RRab (open circles), RRc 
(triangles), RRd (squares), and SPC (filled circles) vari- 
ables. A few of the stars shown in the figure have multiple 
observations. 



3 short period Cepheids in the LMC were analyzed as be- 
ing normal horizontal branch stars. Indeed, the surface 
gravity of these supra-HB stars should be similar to that 
of the RR Lyrae variables, because they are brighter but 
also likely more massive. We then expect that the metallic- 
ity calibration we use, based on normal horizontal branch 
stars, should be roughly correct also for them. In order 
to check this assumption in Figure [221 we show the run 
of the < H > values versus intrinsic (B— V)o colours for 
the LMC variables. These colours were derived from the 
B,V light curves (Di Fabrizio et al. 2004), dereddened 
according to E(B - F)^=0.116 and E(B - F) B =0.086 
(Clementini et al. 2003a), and correspond to the colour 
at the phase of each individual spectrum. Only data of 
stars observed around minimum light are displayed in 
Figure|221 Different symbols are used for the various types 
of variables. Two of the RR Lyrae stars (Al-19711 and A2- 
25301) deviate from the general trend, we suspect that the 
ephemerides and/or the reddening of these two variables 
may be wrong. The short period Cepheids (filled symbols 
in Figure E2J follow the general < H > -(B - V)o re- 
lationship defined by the RR Lyrae variables. The SPC 
spectra closer to minimum light fall on the upper enve- 
lope at redder colours of this relation. This suggests that 
the metallicity calibration we have adopted holds also for 
the SPCs. However, star Bl-05952 falls outside the range 
defined by the RR Lyrae variables, in a totally extrap- 
olated region of this relation, thus its metal abundance 
could be slightly overestimated. 

The [Fe/H] value we obtain for the u> Cen candidate 
Anomalous Cepheid #73 is [Fe/H]= — 1.50±0.12, in agree- 
ment with previous estimates (see Table 6). This star ap- 
pears to be marginally more metal-rich than the main 
metal-poor component of cu Cen at [Fe/H]~ —1.6 (Norris 
et al. 1996), but definitely metal-poorer than both the 
secondary metallicity peak at [Fe/H]~ —1.2 (Norris et al. 
1996), and the metal-rich peak at [Fe/H]~ —0.6 (Pancino 
et al. 2002). 
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For the 3 short period Cepheids in our LMC 
sample we derive: [Fe/H]=-1.96 ± 0.16 for A2-09604, 
[Fe/H]=-1.66± 0.09 for A2-10320, and [Fe/H]=-1.59 : 
±0.03 for Bl-05952, and a corresponding average value of 
[Fe/H]=— 1.74 ± 0.11. Both individual and average metal- 
licities of these three variables are well below the average 
metal abundance of the RR Lyrae stars in our sample 
([Fe/H]=-1.48 ± 0.03) thus suggesting that these three 
short period Cepheids are more likely to be Anomalous 
Cepheids with masses M ~ 1.3M© rather than the short 
period tail of the LMC Classical Cepheids. 

It would be important to derive metallicities for more 
Anomalous Cepheids in different environments. More ex- 
tensive samples of Anomalous Cepheids with adequate 
high resolution spectroscopic data would be clearly wel- 
comed. 

Acknowledgements. A special thanks goes to P. Montegriffo 
for the use of his software, to A. Lay den for providing us up- 
dated ephemerides for the NGC 3201 variables in advance of 
publication, and to A. Piersimoni for sending us finding charts 
for the NGC 3201 variables. We warmly thank the referee A. 
Walker for his constructive and helpful comments. This re- 
search has made use of the SIMBAD data base, operated at 
CDS, Strasbourg, France; it was funded by COFIN 2002028935 
by Ministero Universita e Ricerca Scientifica, Italy 

References 

Alcock, C. et al. 1996, AJ, 111, 1146 

Blazhko, S. 1907, Astron. Nachr., 175, 325 

Bono, G., Caputo, F., Santolamazza, P., Cassisi, S., & 

Piersimoni, A. 1997, AJ, 113, 2209 
Bragaglia, A., Gratton, R.G., Carretta, E., Clementini, G., Di 

Fabrizio, L., & Marconi, M. 2001, AJ, 122, 207 
Brocato, E., Castellani, V., & Ripepi, V. 1994, AJ, 107, 622 
Butler, D., Dickens, R.J., & Epps, E. 1978, ApJ, 225, 148 
Cacciari, C. 1999, in Harmonizing Cosmic Distance Scales in 

a Post-Hipparcos Era, ed. D. Egret & A. Heck, ASP Conf. 

ser. 167, 140 

Cacciari, C, & Clementini, G. 2003, in Stellar Candles for 
the Extragalctic Distance Scale, Lecture Notes in Physics 
(Berlin: Springer), eds. D. Alloin & W. Gieren, Vol. 635, p. 
105 (astro-ph/0301550) 

Carretta, E., Gratton, R.G., Clementini, G., & Fusi Pecci, F. 
2000, ApJ, 533, 215 

Caputo, F., Castellani, V., Marconi, M., & Ripepi, V. 2000, 
MNRAS, 316, 819 

Clement, M.C., Ferance, S., & Simon, N.R. 1993, ApJ, 412, 
183 

Clementini, G., Gratton, R., Bragaglia, A., Carretta, E., Di 
Fabrizio, L., & Maio, M. 2003a, AJ, 125, 1309 

Clementini, G., Held, E.V., Baldacci, L., & Rizzi, L. 2003b, 
ApJ, 588, 85 

Di Fabrizio, L., Clementini, G., Maio, M., Bragaglia, A., 

Carretta, E., Gratton, R.G., Montegriffo, P., & Zoccali, 

M. 2004, A&A, in preparation 
Dolphin, A.E., Saha, A., Claver, J., Skillman, E.D., Cole, A. A., 

Gallagher, J.S., Tolstoy, E., Dohm-Palmer, R.C., & Mateo, 

M. 2002, AJ, 123, 3154 
Fernley, J. A, Carney, B.W., Skillen, I., Cacciari, C, & Janes, 

K. 1998, MNRAS, 293, L61 



Gallart, C, Freedman, W.L., Aparicio, A., Bertelli, G., & 

Chiosi, C. 1999, AJ, 118, 2245 
Gratton, R.G. 1989, Rome Astron. Obs. Preprint 29 
Gratton, R.G., Tornambe, A., & Ortolani, S., 1986, A&A, 169, 

111 

Harris, W.E. 1996, AJ, 112, 1487 

Jurcsik, J. 1995, AcA, 45, 653 

Jurcsik, J. & Kovacs, K. 1996, A&A, 312, 111 

Kaluzny, J., Kubiak, M., Szymahski, A., Udalski, A., 

Krzemihski, W., & Mateo, C. 1997, A&AS, 125, 343 
Kovacs, K., & Walker, A.R, 2001, A&A, 371, 579 
Marconi, M., Fiorentino, G., & Caputo, F. 2004, A&A, in press 

l |astro-ph/0401332| | 
Nemec, J.M., Nemec, A.F.L., & Lutz, T.E. 1994, AJ, 108, 222 
Norris, J.E., Freeman, K.C., & Mighell, K.L. 1996, ApJ, 462, 

241 

Pancino, E., Pasquini, L., Hill, V., Ferraro, F.R., & Bellazzini, 

M. 2002, ApJ, 568, 101 
Pritzl, B. J., Armandroff, T.E., Jacoby, G.H., & Da Costa, G.S. 

2002, AJ, 124, 1464 
Preston, G.W. 1959, ApJ, 130, 507 

Rey, S.-C, Lee, Y.-W., Joo, J.-M., Walker, A., & Scott, B. 2000, 
AJ, 119, 1824 

Rich, M., Corsi, C.E., Bellazzini, M., Federici, L., Cacciari, C. 

& Fusi Pecci, F. 2001, in IAU Symp. 207, Extragalactic 

Star Clusters, ed. E.K. Grebel, D. Geisler, & D. Minniti 

(San Francisco: ASP), p. 140 
Rosenberg, A., Piotto, G.P., Saviane, I., & Aparicio, A. 2000, 

A&AS, 144, 5 
Sandage, A. 1981a, ApJ, 224, L23 
Sandage, A. 1981b, ApJ, 248, 161 
Sandage, A. 1990, ApJ, 350, 603 
Sandage, A. 1993a, AJ, 106, 703 
Sandage, A. 1993b, AJ, 106, 719 

Sawyer Hogg, H. 1973, Publ. David Dunlap Obs. 3, No. 6 
Smith, H.A., Silbermann, N.A., Baird, S.R., & Graham, J. A. 

1992, AJ, 104 1430 
Sturch, C. 1966, ApJ, 143, 774 
Suntzeff, N.B., & Kraft, R.P. 1996, AJ, 111, 1913 
Taribello, E. 2003, Master of Science Thesis, University of 

Bologna 

Udalski, A., Kubiak, M., & Szymanski, M. 1997, AcA 47, 319 

Walker, A.R. 1992a, ApJL, 390, L81 

Walker, A.R. 1992b, PASP, 104, 1063 

Walker, A.R. 1994, AJ, 108, 555 

Walker, A.R. 1998, AJ, 116, 220 

Zinn, R., & West, M.J. 1984, ApJS, 55, 45 



This figure "Al.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



This figure "Bl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



This figure "A2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



This figure "B2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



This figure "A3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



This figure "B3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



This figure "A4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



This figure "B4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405412vl 



